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ABSTRACT 

Deep Chandra observations of the Hydra A Cluster reveal a feature in the X-ray surface brightness that 
surrounds the 330 MHz radio lobes of the AGN at the cluster center. Surface brightness profiles of this feature 
and its close association with the radio lobes argue strongly that it is a shock front driven by the expanding 
radio lobes. The Chandra image also reveals other new structure on smaller scales that is associated with the 
radio source, including a large cavity and filament. The shock front extends 200 - 300 kpc from the AGN at 
the cluster center and its strength varies along the front, with Mach numbers in the range ^ 1.2- 1.4. It is 
stronger where it is more distant from the cluster center, as expected for a shock driven by expanding radio 
lobes. Simple modeling gives an age for the shock front ~ 1.4 x 10** y and a total energy driving it of ~ 10^' 
erg. The mean mechanical power driving the shock is comparable to quasar luminosities, well in excess of 
that needed to regulate the cooling core in Hydra A. This suggests that the feedback regulating cooling cores is 
inefficient, in that the bulk of the energy is deposited beyond the cooling core. In that case, a significant part of 
cluster "preheating" is a byproduct of the regulation of cooling cores. 

Subject headings: cooling flows - galaxies: clusters: individual (Hydra A) - intergalactic medium - X-rays: 
galaxies: clusters 



1. INTRODUCTION 

The X-ray emitting gas at the centers of many cooling flow 
clusters is cooler than the surrounding gas and has cooling 
times in the range 10** - lO** y. Despite this, little of t he gas 
' cools below ^ 1 keV (e.g . Peterson et al. 2001; Tamur a et al.l 

; 1200 U iPeterson et al.ll2()0 3: Kaastr a et al 2004). The main 

issue for cooling flows has become to determine what heat 
source makes up for radiative losses from the gas. Many pos- 
sibilities have been proposed, including thermal conduction 
(e.g.|Naravan & Medvedev 2001), energy released by merg- 
ers (e.g. Motl et aL2004 , ) and heating by an active galactic nu- 
; cleus (AGN: lTabor&BinrievlfT993[ lfucker& David 1997I). 
However, maintaining gas for many cooling times, while pre- 
venting it from cooling to low temperatures or being heated 
until its cooling time is comparable to its age, requires a fine 
balance between heating and cooling. Since cooling flows 
are common (|^bian 199-11), they cannot be a transient phe- 
nomenon, and plausible mechanisms for maintaining the gas 
with short cooling times must almost certainly in volve feed- 
back. While they face some difficulties (e.g. Fabi an et"an 
1200 Ih . AGN heating models provide a natural feedback mech- 
anism, making them the strongest candidate to solve this cool- 
ing flow problem. 

There is clear evidence for some AGN heating in cool- 
ing flows. Large galaxies are prevalent at cluster centers, 
and IB urn si Q390) has shown that cD galaxies at the cen- 
ters of coohng flows have a high incidence of radio activ- 
ity. High resolution X-ray images reveal that radio lobes fed 
by AGN have created cavities in the ho t gas in a growing 
number of clusters (Bohrineer et al. 1993; Carilli et al. 1994; 
i McNamara et al. 2()00: Fabian et al. 200(j: McN amara et al. 
1200 U iBlanton et alJl200U iSchindler et al.ll200U [Heinz et al. 
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120021 lYoung et all l2002h and also in the ho t intersteUar 
medium of individual elh ptical galaxies j Finoguen ov & Jone3 
t200k .Ki-aft et al...200'4l) . IChurazov etah (.2002.) pointed out 
that the enthalpy of a cavity (bubble) can be thermalized 
as it rises buoyantly, providing a heat source for the gas. 
Simple arguments show that the enthalpy is thermalized in 
the wake of a rising bubble, provided that the bubble is not 
too large compare d to its distance from the cluster center 
jBirzan e t alj|2004j). Using optimistic estimates for the heat- 
ing rate.^Bi rzanetal] J2004 ) analyzed data for all cavities they 
could identify in the Chandra archive and found that, while 
the heating rate due to bubbles can be significant, it typically 
falls several times short of radiative losses. 

Standard mode ls of jet-fed radio lobes (e.g. IScheuedll 974t 
iHeinz et al.ll998l i) posit the existence of a "cocoon shock" that 
surrounds the rapidly expanding radio lobes. Although there 
is commonly a bright rim around the radio lobe cavities in 
X-ray clusters, this is composed of coo l (low entropy) gas , 
close to local pressure equilibrium (e.g. iNulsen et al.ll2002l) . 
so that most cavities are not strongly over pressured (Cen A 
is a noteworthy exception, Kraft etal. 2003). However, evi- 
dence has begun to emerge for weak shocks generated by the 
radio lobes in cooling flows. Fabian et al. (2003) found rip- 
ples in the surface brightness of the Perseus cluster, which 
they interpret as weak shocks (sound waves) generated as the 
expansion rate of the lobes varies. They argue that heating 
due to viscous dissipation of sound can make up for radiative 
losses from the Perseus cooling flow. This process has been 
simulated by Ruszkowski et alj (l2004l) . Young et al. (200^ 
found a weak shock surrounding M87 in the Virgo Cluster, 
and Forman et al. (2004) found evidence of further shocks in 
a deeper X-ray image. They determined the energy of the 
outburst that drove the best defined shock, using the result to 
show that heating due to weak shocks is sufficient to stop the 
cooling flow in M87. 

Here we report the discovery of a weak shock front gen- 
erated by an AGN outburst in the Hydra A Cluster Hydra 
A has well known radio lobe cavities (McNamara et al. 2000t 
iDavid et al.l200lHNulsen et al.l2002h . which surround die 1.4 
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GHz inner radio lobes and extend to ^ 40" (40 kpc) from the 
cluster center. By contrast, the weak shock front discussed 
here extends up to ^ 6' from the cluster center, surrounding 
the much larger low frequency radio lobes (Lane et al. 200^ 
and representing a much more energetic phenomenon. 

In section|2lwe give details of the observations and data re- 
duction. In section|3]we discuss the main features of the new 
deep image of Hydra A and in section 0] we discuss the evi- 
dence for the shock front and estimate the energy of the out- 
burst that drove it. The implications of this shock for cooling 
flows and clusters in general are discussed in section |5] We 
assume flat ACDM, with Hq = 70 km s"' Mpc"' and fl^ = 0.3, 
giving a scale of 1 .05 kpc arcsec"' at the redshift of Hydra A, 
z = 0.0538. 

2. OBSERVATIONS AND DATA REDUCTION 

Hydra A was observed with Chandra for 98.2 ksec on 13 
Jan 2004, with ACIS-S at the aim point in VFAINT mode 
(OBSID 4969). A z-sim offset of 5 mm was used to keep the 
cluster center close to the optical axis, while extending the 
field of view to the north to c over the northern low frequency 
radio lobe ( Lan e et al ] l200l . The event list was screened to 
removed ASCA grades 1, 5 and 7, and bad pixels in the stan- 
dard manner. A large background flare made it necessary 
to set the mean background count rate manually when us- 
ing lc_clean to remove periods of high particle background, 
but the mean count rate in ACIS SI after cleaning is close 
to the expected value^. After cleaning, 66.4 ksec of good 
data were left. The data were processed to correct for time 
dependence of the ACIS gain^. Data were filtered accord- 
ing to the prescription of Vikhlinin^ to reduce particle back- 
ground. Background event files were created by processing 
the standard ACIS background files in the same manner as 
the data. Point sources were identified manually for removal 
from spectra and surface brightness profiles. When extracting 
spectra, ARF's and RMF's for extended regions are weighted 
by the number of events in each subregion. ARF's are cor- 
rected to allow for the reduction in low energy response due 
to the build up of contaminant on the ACIS filters. 

3. X-RAY STRUCTURE OF HYDRA A 

Fig- n shows a background-subtracted, exposure-corrected 
and smoothed image of Hydra A in the 0.5 - 7.5 keV band. In 
order to make structures visible over a wide range of radius, 
the image has been divided by a spherical beta model, with 
a core radius of 98", (3 = 0.6 and centered on the AGN. The 
field is covered by ACIS S3 and part of S2 (to the northeast). 
For scale, the chips are 8'.4 wide. The southwest cavity of 
iMcNamara et al. ( 2000) can be seen ~ 0'.5 to the southwest of 
the bright AGN. The corresponding radio cavity to the north- 
east is also visible, though less clearly defined. A 1' bright 
filament stretches from the inner cavity to about the center of 
a larger outer cavity in the northeast (cavity C). The associ- 
ation between the 1' bright filament to the northeast and the 
radio jet of Hydra A is reminiscent of the southwest filament 
in M87 JForman et al.»200 4). suggesting a similar origin. The 
filament in Hydra A is roughly twice as long as that in M87, 
but the radio source in Hydra A is a lot larger too. Cavity C 
has a radius ~ 0'.6 (40 kpc) and occurs at the sharp bend in 
the radio jet leading to the northern 330 MHz radio lobe (see 

^ http://cxc.harvard.edu/contrib/maxim/acisbg/data/README 
* http://hea-www.harvard.edu/~ alexey/acis/tgain/ 
^ http://cxc.harvard.edu/cal/Acis/Cal_prods/vfbkgmd/ 



Fig. 1 . — Chandra 0.5 - 7.5 keV image of Hydra A. The image from the S2 
and S3 chips has been background-subtracted, exposure-corrected, smoothed 
with a 2" gaussian and divided by a spherical beta model centered on the 
AGN. The chips are 8'. 4 wide. The cavities of McNamara et al. ( 2000) can 
be seen close to the blight AGN. The cavity 1' - 2' northeast of the AGN 
(cavity C), at the end of the bright filament, occurs where the 330 MHz radio 
jet bends sharply to the northwest (see Fig.|2). The shock front surrounding 
the radio lobes of Fig.|2]is visible, particularly to the east and north. 
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Fig. 2. — X-ray image of Hydra A with 330 MHz radio contours. The 
image of Fig.Qis shown with logarithmically spaced contours from the 330 
MHz radio map. The two closed contours near the center of the map are at 
the position of the 1 .4 GHz radio lobes, marking the locations of the inner 
X-ray cavities that can be seen in Fig.fTI 



Fig.|3- There is also a deficit in the X-ray emission ~ 2' to 
the southwest of the AGN, associated with the southern 330 
MHz radio lobe. 

There is an edge in X-ray surface brightness, running from 
^ 4'.3 east to ~ 6' north of the AGN, that continues around 
the AGN. The norther n part of t his edge was noted previ ously 
bv lForman et alJ (1200 0) and bv iMarkevitch et al.1 (I2002^ . Al- 
though it is less distinct to the west, the feature can be traced 
all the way around the AGN, except where it falls off the de- 
tectors in the south. We will interpret the feature as a shock 
front and refer to it this way from here on. Fig. |2] shows the 
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Fig. 3. — Surface brightness profiles of tlie sliock front in Hydra A. Upper 
panel: Tlie radial surface brightness profile measured in the sector from PA 
240° to 300° , to the west of the AGN, and in the energy range 0.6-7.5 keV. 
Surface brightness errors are 1 sigma statistical errors. Radial error bars show 
the limits of the bins. The smooth curves show surface brightness profiles for 
shock models with Mach numbers of 1.15, 1.19 and 1.23, from bottom to 
top. Models are scaled to match the observed surface brightness outside the 
shock. Lower panel: 0.6 - 7.5 keV surface brightness profile of the shock 
front in the northeast, measured in a sector from PA 20° to 50° and centered 
2'.9 north of the AGN. Models are shown for shocks with Mach numbers of 
1.26, 1.34 and 1.42. 



same X-ray image as Fig.[T]together with contours of the 330 
MHz radio map of lLane et al The northern radio lobe 

Hes close inside the northern part of the shock front and has 
a similar shape. Although the X-ray image of the shock front 
is truncated to the south, enough of it is visible to see that, 
along with the 330 MHz radio lobe, it lies closer to the AGN 
in the south. The correspondence between the shapes of the 
radio lobes and the shock front supports the interpretation of 
this surface brightness feature as a shock front driven by the 
outburst that formed the radio lobes. 

4. SHOCK MODELS FOR HYDRA A 

4.1. Shock brightness profiles 

The upper panel of Fig.|3]shows the surface brightness pro- 
file of the region to the west of Hydra A, in a sector cen- 
tered on the AGN, from PA 240° to 300°. Point sources 
were eliminated, background subtracted and the resulting pro- 
file exposure corrected. The abrupt change in slope at a ra- 
dius of 202" (211 kpc) is due to the feature that we identify 



as a shock. Outside this, the surface brightness is well fit- 
ted by the power law, r"" with a = 2.65 ±0.17 (90% con- 
fidence). A surface brightness profile was also extracted in 
the region to the northeast of the AGN, where the edge in 
the surface brightness is most prominent. This is shown in 
the lower panel of Fig. |3] In order to match the curvature 
of the radial bins to that of the front, the sector was cen- 
tered at a point approximately 2'.9 north of the AGN (at J2000 
RA, (5 = 9'^18'"07^.30, -12°02'48".4), with the range of PA 20° 
to 70°. Although the data are noisy, the edge is evident in the 
surface brightness profile at a radius of about 200" (measured 
from the center of the sector). 

The break in surface brightness in the northeast is greater 
than that in the west, showing that the shock is stronger in 
the northeast. This is consistent with expectations for radio 
lobe driven shocks. Shock strength is determined by the ratio 
of postshock to preshock pressure. The preshock pressure is 
significantly lower in the north, where the shock is farthest 
from the cluster center. To a first approximation, the pressure 
rise behind the shock is uniform around the front (assuming 
that jet momentum is insignificant on such large scales). Thus, 
the pressure jump will be greater at points on the shock front 
that are further from the cluster center. 

The shock front lies 3'.4 west of the AGN and 4'.3 to 
the east, implying a significant east-west asymmetry in shock 
speed. This is most likely due to asymmetry in the medium 
carrying the shock (as opposed to the outburst driving it) and 
probably the result of higher gas density to the west. If the 
cluster were spherically symmetric, the shock front would be 
running through denser gas in the west than in the east. Any 
pre-existing density asymmetry that slowed the propagation 
of the shock to the west would augment this density differ- 
ence. If postshock pressure is much the same in the east and 
west, then the shock will be weaker in the west, where it prop- 
agates through higher pressure gas. The western shock also 
appears superimposed on a higher "background" of emission 
from undisturbed cluster gas. Together, these effects make the 
shock less distinct in the X-ray image to the west than to the 
east, but it is clearly visible in the upper panel of Fig.|3l 

4.2. Hydrodynamic models 

To quantify the shock, we have used a spherically symmet- 
ric, hydrodynamic model of a point explosion at the center of 
an initially isothermal, hydrostatic atmosphere. For the shock 
front in the west, the initial gas density profile is assumed 
to be a power law, p{r) oc r"**, with rj = 1.82, which makes 
the surface brightness profile of the undisturbed gas consis- 
tent with the observed surface brightness profile beyond the 
shock (see above). The gravitational field (g oc 1 /r) and gas 
temperature are scaled to make the undisturbed atmosphere 
hydrostatic. Surface brightness profiles are determined from 
the model, assuming that the temperature of the unshocked 
gas is 4 keV. The Chandra 0.6 - 7.5 keV response was com- 
puted using XSPEC, with detector response files appropriate 
for these observations and an absorbed mekal model. The 
foreground column density was set to 4.94 x 10^" cm"^, the 
redshift to 0.0538 and the abundance to 0.3 times solar, as ap- 
propriate for Hydra A (results are insensitive to these parame- 
ters, including the preshock temperature). The shock weakens 
as the hydrodynamic model evolves and, since the initial con- 
ditions are self-similar, the flow can be scaled radially to place 
the model shock at the location of the observed shock. Sur- 
face brightness is scaled to match the observed profile in the 
unshocked region. 



4 



The three lines in the upper panel of Fig. |3l show model 
surface brightness profiles for shocks with Mach numbers of 
1.15, 1.19 and 1.23 (increasing upward). A Mach 1.2 shock 
gives a reasonably good fit to the observed shock profile to 
the west. Models were also constructed for the shock to the 
northwest. Here the preshock density profile is flatter, requir- 
ing an initial density power law 77 = 1 .37 (largely because the 
shock does not propagate radially in the cluster). Model sur- 
face brightness profiles are shown in the lower panel of Fig.|3 
for Mach numbers of 1.26, 1.34 and 1.42. The Mach number 
of the shock to the northeast ^1.3-1.4, again confirming 
that the shock is stronger in the northeast than in the west. 

Our simple hydrodynamic model is not accurate well be- 
hind the shock. The initial density profile is only well ap- 
proximated as a power law locally. The shock front is clearly 
aspherical. The outer radio lobes lie close behind the shock 
front, so that they still have an influence in driving the shock 
(if not, they would have been left behind the rapidly moving 
front), violating our assumption that the shock front is driven 
by a point explosion. However, since the lobes are not close to 
the shock front everywhere, it would also be necessary to drop 
our assumption of spherical symmetry to make significantly 
more accurate models. Because of this, the added complexity 
of models with energy fed continuously into the center of the 
flow is not warranted. Despite these shortcomings, the mod- 
els do provide reasonable fits to the surface brightness in the 
region of the shock front. In particular, they give reasonably 
accurate measures of the Mach number (the main source of 
uncertainty is the assumption, implicit in the spherical model, 
that we know the curvature of the front along the line of sight). 

4.3. Physical parameters of the shock 

We base our estimates of outburst age and energy on the 
model for the western shock, since the assumptions of spher- 
ical symmetry and a point explosion are more appropriate in 
this region. The age and energy of the model outburst de- 
pend on the temperature and density of the unshocked gas, 
which we determine by deprojection. The sector to the west 
of the AGN, in which the surface brightness profile of the up- 
per panel of Fig. |3lwas measured, was divided into annular 
regions having at least 4000 source counts in the 0.6 - 7.5 keV 
band. One boundary was fixed at the radius of the shock and 
others fitted around this in order to avoid smoothing across 
the shock. Spectra extracted from these regions were used for 
the deprojection. The method is essentially as described in 
Ipavid et al. ( 2001 ). The spectral model for each ring (annular 
region) includes components for all outer shells as well as the 
current one. Only model parameters for the current shell are 
varied when fitting a ring, so that parameters for each spheri- 
cal shell are determined in turn, working inward. The model 
component for each outer shell is weighted by the ratio of the 
volumes of the intersections between that shell, and the cur- 
rent and outer rings. However, weights for the outermost shell 
were determined on the assumption that the gas density fol- 
lows the same power law as the undisturbed gas in the shock 
model (p cx r"' *^^), from the inside edge of that shell to infin- 
ity, so allowing for the remainder of the cluster. Gas density is 
determined from the XSPEC norm of each model component, 
assuming that the gas density is constant in the corresponding 
shell. For the outermost shell, the gas density is determined 
at the center of the shell, under the assumptions used to deter- 
mine the weights for that shell. 

Even with 4000 counts per spectrum, the errors for the de- 
projected temperatures (~ 1 keV) are too large to reveal the 
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Fig. 4. — Western electron density profile of Hydra A. Deprojected electron 
density versus radius in the sector from PA 240° to 300° . The shock jump 
is at 211 kpc. Density error bars are 90% confidence ranges, but some are 
hidden by the data point. 



temperature jump of the shock. For a shock Mach number 
of 1.2, the temperature jumps by ^ 20% at the shock, but it 
falls rapidly behind the shock due to adiabatic expansion. As 
a result, the temperature rise in the shell immediately inside 
the shock is only expected to be ^ 10%. The preshock tem- 
perature determined from the deprojection is approximately 4 
keV and we use this value for further calculations. 

The deprojected electron density profile is shown in Fig.|3 
Electron densities are determined much more accurately than 
temperatures, chiefly because the overall Chandra count rate 
is insensitive to gas temperature in the range of interest. The 
density jump due to the shock at a radius of 21 1 kpc can be 
seen clearly. It is consistent with the ^ 30% jump expected 
for a Mach 1.2 shock. 

Using the deprojected temperature and density of the un- 
shocked gas to normalize the physical parameters of the shock 
model, we find that the age of the outburst is fs = 1.4 x 10^ y 
and its total energy is £5 = 9 x 10^*^ erg. The main source of 
uncertainty in the age of the shock is the preshock tempera- 
ture, since the Mach number is not sensitive to the model (if 
the shocked cocoon is axially symmetric, the size we measure 
is not affected by projection). Shock energy is proportional 
to the preshock temperature. It also depends on the density 
profile of the unshocked gas. There is a substantial geometric 
uncertainty, since the shock front is clearly aspherical. Our 
estimate uses the piece of the shock that is nearest to the clus- 
ter center, underestimating the total shocked volume, so tend- 
ing to underestimate the total energy. Lastly, our estimate is 
affected by the inaccurate assumption of a point explosion. 
Despite all this, we expect the true energy of the outburst to 
be no more than a factor ^ 2 astray. The mean mechanical 
power of the outburst is = E^/t^ ~ 2 x lO'*^ erg s"'. 

4.4. Shock vs Cold Front 

Interpreting the shock feature as a cold front instead would 
be difficult, although temperatures are not determined accu- 
rately enough from the current data to categorically rule this 
out. The density discontinuity (Fig. |4} shows that there is a 
front. If it is a cold front, the pressure must be continuous 
there. Otherwise the only physically reasonable interpreta- 
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Fig. 5. — Western temperature profile of Hydra A. Projected temperature 
versus radius in the section from PA 240° to 300° . Error bai's shows 90% 
confidence ranges for the temperature. 



tion is a shock. For a density decrease ^ 30%, the tempera- 
ture must increase by ~ 30% to keep the pressure continuous. 
Fig-HJshows projected temperature versus radius for the sec- 
tor where the deprojection was done. Temperatures were ob- 
tained by fitting absorbed mekal models, with the absorbing 
column fixed at the foreground value (A^h = 4.94 x 10^" cm~^) 
and the abundance free. If anything, temperature decreases 
with radius at the front (r = 2 1 1 kpc). The same is true at other 
PA's, where greater density steps would require greater tem- 
perature steps. According to the shock model, the maximum 
rise in projected temperature should be ~ 5%, confined to the 
two rings behind the shock. This is consistent with Fig. |5l 
although there is a mild negative temperature gradient in the 
vicinity of the shock that is not included in our simple shock 
model. The association between the front and the radio lobes 
is much better explained by a shock driven by an outburst 
than by a cold front. Furthermore, cold fronts are thought to 
be contact discontinuities formed when a body of gas moves 
through gas of higher entropy. In that case, it would be diffi- 
cult to form a closed cold front, whereas a shock driven by an 
outburst at the cluster center is expected to surround the AGN. 
Altogether, the evidence strongly favors interpretation of this 
feature as a shock front. 

5. DISCUSSION 
5.1. Jet History 

X-ray measurements of shock properties provide a new 
tool for the study of AGN outbursts and radio sources. As 
discussed above, the uncertainty in the age of the outburst, 
~ 15%, is mainly due to uncertainty in the the preshock tem- 
perature. Shock strength is determined by the overpressure 
of the radio lobes. If it remains constant, then, under the 
assumptions of our simple shock model, the total energy of 
the shocked region scales as £5 k rl~^ oc fs^ '', where rj is the 
shock radius and fs is the time since the outburst. Jet power 
needs to vary with time as t'^~^ to keep up with this. The shock 
weakens or strengthens with time, depending on whether jet 
power lags or leads this time dependence. The fraction of the 
total energy going into cavity enthalpy also depends on the 
overpressure of the cavities. Detailed modeling will constrain 



the history of the AGN outburst. 

Although there are two X-ray cavities visible along the 
northeastern radio jet, none is evident at the location of the 
northern radio lobe. However, the surface brightness is low 
there and has structure on the scale of the lobe (including the 
shock front) that helps to mask any cavity if it i s present. As 
found for the inner cavities (Nulse n et alJl2002l) . cavity C in 
the northeast is made easier to recognize by the bright rim 
along its outer edge (Fig.[0. Without this, the contrast of cav- 
ity C would only be ^ 10%, making it considerably harder to 
see. 

The structure of the northeastern jet and the shock front 
suggest a complex history for the radio source. If the bulk 
of the northeastern jet was filled with a connected body of 
radio plasma at the onset of the current outburst, a compres- 
sion wave (shock) would have been launched simultaneously 
from its whole length. If, as it is now, the shock was strongest 
in the regions of lowest preshock pressure, near to the outer 
lobes, it would have propagated farthest in th ose regions (as 
in the simulations of Ruszkowski et al. 2004). This is incon- 
sistent with the bulges in the shock front to the east and west 
of the center, suggesting that this did not happen. At the other 
extreme, the outburst could have launched a new jet that has 
penetrated ~ 300 kpc from the AGN to the northern radio lobe 
in 1.4 X 10*^ y, with a mean speed of vj ~ 2100 km s"'. It 
appears unlikely that jet momentum would have been signifi- 
cant for the propagation of the jet beyond the location of cav- 
ity C, where it turns through ~ 90°. If so, the mean speed of 
the lobe is excessive f or a bubble that was dr iven most of the 
way by buoyancy (e.g.'Churazov et al."200r). The intermedi- 
ate case, of a newly launched jet connecting after a time into 
aging trails of radio plasma left by earlier outbursts can avoid 
these two extremes. It also requires more modest asymme- 
tries in gas properties than the pure jet model to account for 
the north-south asymmetry of the radio structure and shock 
front. 

5.2. Cooling Flow Feedback and Preheating 

It needs to be asked whether the outburst in Hydra A has 
any connection to the cooling flow. If so, then the gas that 
fuels it must be cooled or cooling intracluster medium (ICM). 
The line emitting gas at the center of Hydra A might have 
been deposited by a merger (iMcNamarall 1 995h . which could 
also be the source of fuel for the outburst. As discussed in 
the introduction, the strongest argument for a connection be- 
tween the outburst and the cooling flow is the difficulty of 
maintaining cool cores without feedback, but this does not 
rule out occasional mergers leading to larger than usual out- 
bursts. On the other hand, large outbursts appear to play 
a significant role in balancing radiative losses (see below). 
In that case, it would be surprising if they were not part 
of the feedback process. The correlation betwe en cooling 
luminosity and cavity heating power found by Birzan et alJ 
( 2004) also suggests strongly that the two are connected by 
feedback across the full range of systems. The outburst in 
Hydra A is three orders of magnitude more energetic than 
that which produced the 14 kpc shock in M87 (^ lO*"*^ erg, 
iForman eTa l. 2004). The outbursts responsible for the ripples 
in Perseus (Fabian et al. 2003) lie somewhere betw een these 
in energy. MS0735.6H-7421 (McNamara et alJ2004 ) and Her- 
cules A (Nulsen et al., in preparation) both have more ener- 
getic outbursts than Hydra A. At the least, outbursts span a 
large range of energies. A large sample is required to deter- 
mine whether the large outbursts form a distinct class from 
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those required to regulate cooling cores. 

In the absence of heat sources, radiating gas cools to low 
temperature in approximately one cooling time, even after al- 
lowing for flow induced by the cooling. Thus, if the mass of 
gas inside radius r is Mg(r) and the cooling time of gas at ra- 
dius r is tc(r), a good approximation to the cooling rate when 
the gas between ri and is cooling to low temperatures is 
jFabian&Nulsenll979h 

^■j^_ M,{r2)-M,{n) 
tc(r2)-tciri) 

Applying this to the results of iDavid et al.l (1200 Ih for Hydra 
A shows that the cooling rate would reach 100 M0 y"' al- 
most immediately at the onset of cooling, and climb steadily 
thereafter Such high cooling rates have yet to de monstrated 
convincingly for any cluster (e.g. iPeterson et ani200 3) and 
so must be rare. This requires that, if there is any period 
when the gas cools to low temperatures, it must be brief com- 
pared to observed central cooling times. Thus another out- 
burst must occur before, or v ery soon after, the central gas 
begins to cool in 4 x 10^ y (iDavid et al.ll2001h . The max- 
imum time between the last outburst and the next is there- 
fore fmax 5 X 10^ y, and the minimum average power re- 
quired for such outbursts to prevent the gas from cooling is 
Pniin — £s/fmin — 6 X 10"^^ erg s"', significantly more than the 
2.5 X 10** erg s"' radi ated fr om within 150 kpc of the clus- 
ter center fDavid et aTl l200l . The outburst in Hydra A can- 
not be a very efficient process for regulating the cool core. 
The shock front and outer radio lobe already lie outside the 
cooling flow region, "wasting" energy in the cluster at large 
radii and making the inefficiency manifest. Energy that is de- 
posited outside the cooling flow r egion adds to the t otal en- 
ergy of the ICM and, as noted by David et al. ("2001'), it can 
make an appreciable contribution to cluster "preheating." For 
a cluster temperature of 4 keV, the virial mass of Hydra A 
is about 5.4 x lO'"* Mq (iBrvan & Normanlll998h and its gas 
mass about 7.6 x 10'-^ Mq (for a gas fraction of 14%). The 
mean heating power required to add 1 keV per particle to this 
gas over lO"* y is 8 x 10'*"^ erg s"'. Our estimates of the mean 
heating power from AGN outbursts lie in the range 6 x 10"''* 
- 2 X lO"*^ erg s~', so that the "waste" power could easily be 
a major source of preheating in Hydra A. There is no rea- 
son to believe that outbursts in much smaller systems (e.g. 
IFinoguenov & Jones 2001) are more efficient, so these may 
also make significant contributions to preheating. 

■Birzan et al. found that bubble heating alone is in- 

sufficient to make up for radiative losses in most cooling 
flows. It should be noted that Hydra A was one of the ex- 
ceptions. The total of pV for t he two 1 .4 GHz i nner cavities 
in Hydra A is 1 .7 x lO^** erg (Birzan et al. 2004). Taking the 
pressure from David et al. (J^OOl), for cavity C to the north- 
east we estimate pV ~ 3.7 x 10^^ erg. Doubling this to allow 
for outer cavities to the southwest, the total of pV for cav- 
ities in Hydra A is now ^ 9 x 10^'' erg. If the cavities are 
filled with relativistic plasma, their total enthalpy is the sum 
of 4pV for each, roughly one third of the estimated total en- 
ergy of the shock. The shocks in Hydra A and other systems 
dFabian et alJ2003llForman et al I200I show that bubble en- 
thalpy may not be the most significant form of mechanical 



energy produced by AGN outbursts. As noted above, two 
systems have been found with even more energetic shocks 
(Hercules A, Nulsen et al., in preparation; MS0735.6H-7421, 
McNamara et al. 2004). If they are a general feature of AGN 
outbursts, shocks will certainl y help to account for the en- 
ergetics of cooling flows. Birzan_et al,! (2004) found a large 
scatter in the ratio of radio power to bubble power, suggesting 
that current radio power is not a good measure of the history 
of AGN outbursts, even on timescales ~ 10^ - 10** y. How- 
ever, while the central galaxies in cooling flow clusters are 
frequently active in the radio, radio sources as large as Hy- 
dra A are exceptional. If such outbursts do occur in a signifi- 
cant proportion of all clusters, they can only be powerful radio 
sources for a relatively small fraction of the time. 

5.3. Black Hole Growth 

The mass, Ms, that was accreted by the AGN to fuel this 
outburst can be related to the shock energy by = eM^c^, 
where e is the efficiency of the AGN for converting accreted 
mass into jet power (or other mechanical output). With our 
value for the shock energy, this gives Ms ~ 5 x lO^e"' M0. 
This exceeds the esti mated mass of the bl ack hole in Hydra A, 
Mh ~ 4 x 10'' Mq ( Sambr unaetalJ20 00). unless the mechan- 
ical efficiency of the AGN e > lO'l The mean mechanical 
power of the outburst, f s — 2 x lO"'^ erg s"', is in the range 
of quasar luminosities, but a modest fraction of the Eddington 
limit for such a massive black hole. It is well in excess of the 
power needed to make up for radiative losses from the cooling 
flow. 

6. CONCLUSION 

A deep Chandra X-ray image of Hydra A has revealed new 
structure associated with the Hydra A radio source. This in- 
cludes a break in the surface brightness 200 - 300 kpc from 
the cluster center that we interpret as a shock front. The 
shock front surrounds the low frequency radio lobes, closely 
in places, showing a significant resemblance, and we interpret 
it as the cocoon shock of the radio source. 

The Mach number of the shock varies in the range ~ 1 .2 
-1.4 around the front. Simple modeling shows that the 
shock outburst commenced 1.4 x 10** y ago and its energy 
^ 9 X 10*" erg. The mean mechanical power of the outburst 
^ 2 X 10*^ erg s"', in the range of quasar luminosities. Such 
large outbursts can easily account for the radiative losses from 
the cooling flow in Hydra A. If they are part of a feedback cy- 
cle that sustains the cool gas in Hydra A, then the feedback 
heating process is inefficient, in the sense that more than half 
of the mechanical energy output of the AGN is deposited in 
regions outside the cooling flow. Over time this will have pro- 
duced significant "preheating" of Hydra A. 
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